Dipole-Dipole interaction in photonic crystal nanocavity 
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Dipole-dipole interaction between two two-level 'atoms' in photonic crystal nanocavity is investi- 
gated based on finite-difference time domain algorithm. This method includes both real and virtual 
photon effects and can be applied for dipoles with different transition frequencies in both weak and 
strong coupling regimes. Numerical validations have been made for dipoles in vacuum and in an 
ideal planar microcavity. For dipoles located in photonic crystal nanocavity, it is found that the 
cooperative decay parameters and the dipole-dipole interaction potential strongly depend on the fol- 
lowing four factors: the atomic position, the atomic transition frequency, the resonance frequency, 
and the cavity quality factor. Properly arranging the positions of the two atoms, we can acquire 
equal value of the cooperative decay parameters and the local coupling strength. Large cooperative 
decay parameters can be achieved when transition frequency is equal to the resonance frequency. 
For transition frequency varying in a domain of the cavity linewidth around the resonance frequency, 
dipole-dipole interaction potential changes continuously from attractive to repulsive case. Larger 
value and sharper change of cooperative parameters and dipole-dipole interaction can be obtained 
for higher quality factor. Our results provide some manipulative approaches for dipole-dipole in- 
teraction with potential application in various fields such as quantum computation and quantum 
information processing based on solid state nanocavity and quantum dot system. 

PACS numbers: 42.50.Ct, 34.20.-b, 37.30.+i 



I. INTRODUCTION 



Since Purcell predicted spontaneous emission rate 
could be changed by electromagnetic environment in 1946 
[H, the effect of electromagnetic field on radiation prop- 
erties has been thoroughly investigated [2hllol|. and clas- 
sified in the category of cavity quantum electrodynamics 
(QED). Characteristics of QED have been greatly investi- 
gated both theoreticall y an d experimentally [h1416T [ , and 
many kinds of devices [171 [l8[ based on this theory have 
been developed. Concepts such as enhanced and inhib- 
ited spontaneous emission 0, S, reversible spontaneous 
emission Q , photon blockade [1J| , micromasers [l9[ , low 
threshold lasers |2jj,|2l|, etc., have become very familiar. 

Dipole-dipole interaction could also be greatly modu- 
lated by the electromagnetic environment. Photon emit- 
ted by one dipole could be absorbed by the other or vice 
versa. The strength of interaction is decided by pho- 
ton emission, transmission and absorption. Many dif- 
ferent kinds of electromagnetic environment can be used 
to control or change these characteristics, such as vac- 
uum [3L I22T - I25] . optical cavities 26-29|, optical lens [30j j . 
dielectric droplet |3lj . photonic materials 32-38J, metal 
surface [39Tl4a ] , metamaterial 0, SBf and so on. For ex- 
ample, optical lens and waveguide have been designed to 
collect the emission photon and transfer it to the other 
dipole. Optical cavity or metal surface can enhance the 
emission or absorption rate roughly by the ratio of qual- 
ity factor Q and the mode voluum V. 



Photonic crystal nanocavity is one of the promising 
platform to investigate dipole-dipole interaction, because 
local coupling strength of dipole and photon can be tai- 
lored and integrated to photonic crystal waveguide is ex- 
tremely convenient. High quality factor Q — 2.5 x 10 6 
and small mode volume V ~ (A/n) 3 have been realized 
for photonic crystal cavity [46[ . Numerical investigations 
show that ultra-high quality factor with Q ~ 10 9 can 
be designed through finely tuning the scatters around 
the cavity with little change of the mode volume [47j |. 
Furthermore, static and ultra fast dynamic control of 
the cavit y frequ ency and the quality factor have been 
3. 
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On the other hand, temperature 



, electric field 54], magnetic field 55|, 56 1 are 
much suitable for fine tuning the energe levels of dipoles 
located at certain position in solid system. 

Recent studies show that this kind of interaction could 
be used to implement quantum entan glem ent preparation 
and quantum information processing [57H64| . coop erative 
radiation 0J^55| , Forster energe transfer|66l. |67|. dipole 
nanolaser [68|, and so on. Furthermore, some novel quan- 
tum phenomena muz! has been found. All of these ap- 
plications and phenomena are related to the cooperative 
decay parameters or dipole-dipole interaction potential. 
According to Eq. (|5c|) . dipole-dipole interaction potential 
can be got through the cooperative decay parameters. 

In the previous theoretical studies, mode-expansion 
method J^jLpl. H^-Hl or Green function method 
aims |74| is often adopted to investigate this cooper- 
ative decay parameters. These two methods work well for 



electromagnetic environment with perfect boundary con- 
dition. Because of the extremely complexity of finding 
eigen-mode, mode-expansion method can be used only 
for simple case such as vacuum or planar cavity. Be- 
sides, the exact analytic Green function is also hard to 
be obtained for complex electromagnetic environments. 
Numerical method is necessary for studying this kind of 
interaction in photonic crystal nanostructure. 

In this paper, we put forward a simple numerical 
method to investigate dipole-dipole interaction in pho- 
tonic crystal nanocavity through finite-difference time 
domain algorithm (FDTD). We calculate the collective 
and individual radiation rates of classical dipoles by di- 
rectly solving Maxwell's equations in real space with a 
free-space boundary condition. By using the result of 
two dipoles radiation rate minus the sum of the two in- 
dividual radiation rates, we get the cooperative decay 
parameters and dipole-dipole interaction potential. A 
similar method has been used for local density of pho- 
tonic states calculation [75l - l77l | . It is direct for structures 
with arbitrary shape and physical quantities such as radi- 
ation power at different frequencies can be got in a single 
simulation run. Numerical results in vacuum and planar 
conductor cavity show that our method works well for 
this kind of investigation. For dipoles located in pho- 
tonic crystal cavity, the strength of dipole-dipole inter- 
action depends heavily on the following four factors: the 
atomic position, the transition frequency, quality factor, 
mode volume and the cavity frequency. Effect of these 
factors on the interaction strength are also shown. 

This paper is organized as follows: In Sec. II, we give 
the model of dipole-dipole interaction with some illumi- 
nations. Numerical method is presented in Sec. III. Vali- 
dation of this method is presented in Sec. IVA for dipoles 
in vacuum and in planar microcavity. In Sec. IVB, ap- 
plying our method to photonic crystal nanocavity, we in- 
vestigate the effect of the atomic position, the transition 
frequency, quality factor and the cavity frequency on the 
dipole-dipole interaction strength. In Sec. V, we give a 
brief conclusion. 



II. MODEL 



2 




B 



FIG. 1: (Color online). Schematic diagram of dipole-dipole 
interaction. Consider two two-level 'atoms' A and B. Atom A 
(B) has two states, the ground state |(Ja} (|<7b)) and the ex- 
cited state \eA) (|es)), with a transition frequency uua (<^b). 
Atom A is in the excited state \e A ) and atom B is in the 
ground state |<?s). They both interact with electromagnetic 
field with eigen frequency co n . Analogous to the process of 
Lamb shift where both virtral and real photon exchange be- 
tween the atom and the quantum electromagnetic field take 
effect, photom exchange between the two atoms contribute to 
dipole-dipole interaction. 

H = H + V, 

H = h^2 oji\ei)(ei \ + Hymnal a n , 

i=A,B n 
V = h J2 Y}9n^W n \g l ){e l \ + c.c]. (1) 

i=A,B n 

where and a n are, the photonic creation and 
annihilation operators, ui n is the frequency of the 
eigen mode of the electromagnetic field, g n {^i) — 
JWj(2e fio; n ) _1 / 2 E„(ri).Ui ( i = A,B ), is the the cou- 
pling coefficient, E„(r) is the eigen mode, is the tran- 
sition dipole moments of atoms i, and £o is the vacuum 
permittivity. In Eq. ([T]) , Hq is the noninteracting Hamil- 
tonian, and V represents the interaction between dipole 
and electromagnetic field. There are three states for the 
system considered: (1) atom A is in the excited state, 
and atom B is in the ground state, without any photon, 
i.e., | a) = \e A , 3b, 0), (2) atom B is in the excited state, 
and atom A is in the ground state, without any photon, 
i.e., \b) = \gA,eB,0), (3) both atom A and atom B are 
in the ground state, with a photon of frequency u;„, i.e., 

\ c n) = \9A,9B, In)- 

The initial state is prepared in \a). Then, the state of 
the system evolves as 



Schematic diagram of dipole-dipole interaction is illus- 
trated in Fig. 1. There are two two-level 'atoms' A and B 
located at r A and yb respectively and they both interact 
with the electromagnetic field with eigen frequency u) n . 
The Atom A (B) has two states, the ground state \g A ) 
(|<7b)) and the excited state \e A ) (|es)), with a transition 
frequency uj a (wb). The Hamiltonian of the system in 
the rotation wave approach reads [78| : 



\m) =a(t)W)+b(t)\b)+Y,Cn(t)\c n ) = U(t)\a). (2) 

n 

where U(t) is the evolution operator and a(t) — 
(a\U(t)\a), b(t) = (b\U(t)\a), c n (t) = (c n \U(t)\a). To 
derive the equation for the atom-field dynamics nonper- 
turbatively, the resolvent operator G(z) = l/(z — H) is 
adopted [79(. The corresponding advanced and retarded 
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propagators are G ± (E) = lirn__ >0 + G(E ± i-q). Then the 
evolution operator can be expressed by 



U(t) 



+ 00 



dw[G-{u) - G + (u)]exp(-iojt)/2m. (3) 



The dynamic properties of this system are governed by 
G^ a (us) and G^ a . The matrix elements of the resolvent 
reads 



(Z ~ hhJ A )G aa = 1 + ^2 V ac G ca , 
n 

(z - hu B )G ba = ^ V bcGca, 
n 

(z - hw n )G ca = V ca G aa + V cb G ba . (4) 

where G aa = (a\G(z)\a), G B a = (b\G(z)\a), G ca = 

(c nk \G(z)\a) and V ac = V* a = (a\V\c), V bc = V c * b = 
(b\V\c)_ 
Eliminating G ca , we have 



G aa (z) 
G ba (z) 



= (z- hw B - W B b(z))/E, 
= W BA (z)/E. 



(5) 



E is given by 



E = [z- hoj A - W A a] [z - fu^ B - Wbb] - WabWba (6) 

The local coupling between atom and the field 
{Waa,W B b ) or the dipole-dipole coupling between atom 
A and atom B (Wab, Wba) are denoted by 



Wij{z,T U Tj) = 



gn(r«)gn(rj) 



(7) 



Z — huJ n 

Clearly, these terms can be written as 

Wjj^TuTi) = fi[A 13 (.,r 1 ,r 3 )T» r " 3( " 2 r " r,) ], 

(8a) 

Tyfar^Tj) = 2-!r^2g*(r i )g n (r j )8(uj - u n ), (8b) 

n 

A^r^r,) = ± P f +QC ' dz Tij{z ' Ti ' rj) . (8c) 

l-K J Q UJ - Z 

To better understand the physics underlying of these 
equations, we make some illumination. We are able to 
self-consistently determine the roots of E — (Eq. ©). 
The real parts are corresponding to the energy levels of 
these two atoms and the imaginary parts are the life- 
time. Different from the usual Markov approximation 
method where Wij = A, B) are independent of uj and 
have been replaced by their approximate value Wu(u)i) 



and Wij((cJi + ujj)/2), this method is non-Markov and 
this self-consistent process can give the rigorous dipole- 
dipole coupling. We can also make use of the Markov 
approximation values Wij — A,B) in Eq. ([6]) and 
get the same results as usual. However, for complex elec- 
tromagnetic environment such as photonic crystal or pho- 
tonic crystal nanocavity, Wij may vary sharply around 
some frequency and non-Markov is necessary. Owing to 
Wij, equation. ([6]) also shows that both the energy lever 
and the lifetime of the states are spliting. Different cou- 
pling Wij makes different energy lever spliting and differ- 
ent lifetime spliting. Dipole blockade, which have been 
widely investigated in quantum information processing 
needs large energy gap. Efficient superradiant emission, 
steady state entanglement preparation and fast Forster 
energy transfer needs large lifetime spliting. Inspired by 
these novel application, accurate values for Wij are ex- 
pected. Furthermore, equation (J7J shows that all of the 
photonic eigen modes contribute to the local coupling 
and dipole-dipole interaction and both real and virtual 
photon effects have been taken into account. 

Insert g n (r l ) = iu> l (2e Q fuj n )- 1 / 2 E n (r i ).u l (i = A,B) 
into Eq. (|8b[) for i ^ j, and define Si.j(u>) = 
TTu>iUJjUiUj I (eohu) (i = A.B), where and are the 
size and the unit vector of transition dipole u^. Then, 
the cooperative decay parameters reads: 

Tij(u},ri,rj) = s it j(ui)^'En(ri).\i i E n (rj).Uj6(ui - ui n ). 

n 

(9) 

For simplicty, we denote Tij(u>) and A.y(w) for 
Tij(ui,Yi,Yj) and Aij(ui, Tj, r^) respectively. 

From equations 2.12a to 2.14a of reference 80], it is 
easy to see that if {E„(r)} compose a complete set of 
eigen mode, {E* (r)} is also a complete set of eigen mode. 
Then 

Tij(u) = Si,iMy^E*(ri).UjE ra (rj).u i (5(a; - ui n ). (10) 

n 

could also be written as: 
Tij(uj) = s lJ (w)J]E„(r l ).u l E;(r J ).u^(w- w„). (11) 

n 

So we have : 



= ^^5I[E n (r i ).u i E;(r J -).u J - + /i.c.]<J(w-a; n ). 

n 

^ (12) 

Once we get the cooperative decay parameters Tij(uj), 
the dipole-dipole interaction potential (w) can be 
achieved through Eq. (I5cl) . 



III. METHOD 

Here we propose a new method to rigorously get the 
cooperative decay parameters based on finite-difference 
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time domain algorithm. We show that dipole-dipole in- 
teraction can be got through calculating the collective 
and individual radiation rates of classical dipoles. We 
begin with Maxwell equations: 



V x E(r,i) 



dB(r,t) 
d1~' 



. . , ,0E(r,i) dP(r,t) 

V x B(r,t) =Moe(r) — J-^+Mcr 



dt 



Of 



V-e(r)E(r,i)=p(r,t), 
V • B(r,i) = 0. 



(13) 



Expand E(r, t) = J2 n «nWEn(r) where E„(r) is the 
same as the eigen modes at the quantum analysis sec- 
tion. If we let P(r,t) = e _J " ot u(r), then a n (t) satisfies 
the following equation: 



2 — iujQt 



dru(r) • E;(r). (14) 



a n {t) + u n a n (t) = -u 
Through Eq. (fI3)l. we get: 

*n(t)= Um ^H^ /dru(r).E;(r). (15) 
If u(r) = u,5(r— rj), the emission power is given by 



pm = -^i e 1 E a * • E n( r ^) i 2 <^ - ( 16 ) 

n z 

In the two dipoles instance, i.e., i — A, B, the total 
power is 

P^(u, ) = ~^El E Si-E^rOI 2 ^-^). 

n i=A.B 

(17) 

If there is only one dipole, i.e., i = A or i = B, then 
P> ) = Jc 2 El "A- E;(r A ) | 2 5(a; - w n ), (18) 

n 

P B K) = Jc 2 El E:(r fl ) | 2 <S("o - u n ). (19) 

n 

Combining Eq. jT7|, Eq. (JTSJl and Eq. (fTU), we difine 
the cooperative emission power as: 

P co (w ) = P AB {u a ) - P A (lo ) - P B (u ), 

= ^J2^-K(rA)u B -E n (r B )+ 

n 

h.c.}5(u - w„). (20) 
Through Eq. (gDJ) and Eq. ([12]). we found that: 



(21) 



where T^ (w) and P,? (w) correspond to cooperative decay 
parameters and cooperative emission power of two unit 
classical dipoles in complex environment. Pq(uj) corre- 
spond to the emission power of a unit classical dipole in 
vacuum. r° (u)) is the local coupling strength for an two 
level atom with transition dipole moment and tran- 
sition frequency of bare atom uii in vacuum. Using the 
result of [1], we get T^(u;) / (ufojf) = u>/ (3ne hc 3 ). 

In order to make the cooperative decay parameters 
more clearly, wc define 



OLi 



2P »' 



(22) 



where ctj is totally decided by the two level atom i. Then 
we get 



(23) 



Equation. (|23[) is the main content of our method. 
t](lu) can be obtained through calculating P c c D (u;) and 
P£(lu) by FDTD method. Utilizing Eq. (JHcJ), we can 



get the dipole-dipole interaction potential Ay-(o>). Fur- 
thermore, we see that our method may be generalized 
to study many dipole interactions through this similar 
procedure. 



IV. NUMERICAL RESULTS AND DISCUSSION 

A. Cooperative decay parameters in vacuum and 
planar microcavity 

Before applying our method to investigate dipole- 
dipole interaction in photonic crystal nanocavity, we 
first numerically validate our approach in the case for 
dipoles in vacuum or planar microcavity where analytic 
expression for the cooperative decay parameters r?-(cj) 
can be easily got through mode-expansion method. In 
this section, for simplicity, we set u>a — = <^o & n d 



OLA 
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/Sireohc 3 . In vacuum, the two 
transition dipole moments are parallel to the x axis and 
the interatomic separation axis is along the z axis. The 
eigen- function of vacuum reads Ek,z(r) = e k e* kr . Insert- 
ing this into Eq. (fT2|) and making use of Xw=i 2 = 
Sij — kikj, we get the analytic cooperative decay param- 
eters Tij(uj): 



3uTo(u;o) r sin(wP/c) cos(uR/c) sin(ujR/c) 



2w L (uR/c) ' (wP/c) 2 (ljR/c) 3 

(24) 

where R = |r< — Tj\ is the distance between the two 
dipoles. If R = 0, for the two dipole located at the 
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same place, Tij(uio) = Fo(o;o) = ufaf /Sireohc 3 which 
is the radiation rate of the dipole in vacuum. We set 
u = ujq = 2irc/X and increase the interatomic separation 
R to simplify our discussion. In Fig. [2j black line is for 
the analytic results (Eq. ([21)1 ) and red boxed dots repre- 
sent the calculated results (Eq. (|2"5)l ). It can be clearly 
seen that r^w) got through our method agrees very well 
with the analytic results. 

1.0 

0.8 - 
0.6 - 



0.4 - 



0.2 



0.0 - 



-0.2 



-analytic result 
calculated result 



■>= 0.2 




-0.4 I ■ 1 ■ 1 ■ 1 ■ 1 

2 4 6 8 

R/Xo 

FIG. 2: (Color online). Comparison between our numerical 
results and the analytic results of cooperative decay parame- 
ters for dipoles located in vacuum. Cooperative decay param- 
eters Tij/Fo as a function of the interatomic separation R/Xo 
when u = ujq = 2nc/Xo; Black line is for the analytic results 
(Eq. I|24p ) and red boxed dots represent the calculated re- 
sults (Eq. ((23j) ) . The orientation of the two transition dipole 
moments are the same and perpendicular to their separation. 



FIG. 3: (Color online). Sketch of setup and the comparison 
between our numerical results and the analytic results of co- 
operative decay parameters for dipoles located in vacuum, (a) 
Schematic diagram of the planar cavity with a top view of the 
computation domain at the right side. The two dipoles locate 
at the center. The orientation of the two transition dipole mo- 
ments are the same and along z (red arrowhead). They are 
arranged along x axis and their separation is R. (b) Coopera- 
tive decay parameters r^/To as a function of the interatomic 
separation R when z\ = Z2 = L/2, L/Xo = L/X = 0.7 ; black 
line is for the analytic results from Eq. (|25p and red boxed 
dots represent the calculated from Eq. (l23l) . 



Similar result is also got when the two point dipoles 
are located in a planar nanocavity. Here, the cooperative 
decay parameters Tij(uj) could also be got analytically. 
The structure is shown in Fig. G^a). The origin is at 
the center for x, y axis and at the bottom conductor for 
z axis. The CPML encloses the domain, as can be seen 
from the top view at the right-hand side. The orientation 
of the two transition dipole moments are the same and 
along z (red arrowhead) . They are arranged along x axis 
and their separation is R. The analytic cooperative decay 
parameter Tij(u>) for these two dipoles is: 



3A 



,2-kR 



Tij[u) =r M(^){ / d0oos(— cos0)+ 



,nX 



X 



[2L/X] 

2[1 - (t^) 2 ] cos(-^) cos(-^)}x 

71=1 ~ " J ' 



J 2 J dQ 008(2^(1)2 COS 



0). (25) 



where ui = 2irc/\, [2L/X] is the largest integer less than 
2L/X, rrj(u;n) is spontaneous emission rate in vacuum for 
any one of the two same dipoles, and Zi is the coordinate 



of z for dipole i, here, we set Z\ = z 2 = L/2, R is the 
distance between the two dipoles in xy plane. For R = 0, 
this leads to the well-known results 



r,( Wo ) = r ( Wo )(^){i- 



[2L/X ] 

E 

71=1 



2[1-(^)> S 2 (^)} 



2L 



77,71" N 

(26) 



In the calculation, we set L/Xq = L/X = 0.7, and 
R/L is gradually increased from to 49/41 for simplicity. 
Cooperative decay parameters T^/To as a function of 
the interatomic separation R is displayed in Fig. E^b). 
black line is for the analytic results from Eq. (|25p and 
red boxed dots represent the calculated from Eq. (j23j) . 
We find that the numerical result from our method also 
agrees very well with the analytic solution. 
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FIG. 4: (Color online). The sketch of the photonic crystal 
nanocavity. It consists of a thin dielectric slab with air holes 
arranged as triangular lattice. The lattice constant is a, the 
radii of the air hole is r — 0.3a and the slab height is d = 0.6a. 
The refractive index of the slab is 3.4. There is a defect hole 
in the center with refractive index ride/ = 2.4. Its radii is 
the same as the air hole. The origin of the axes is set at the 
center of the photonic crystal nanocavity. The two dipoles 
are located at point A and B on the center plane of the slab 
and the two dipole moments are parallel to the x axis. 



B. Dipole-dipole interaction in photonic crystal 
nanocavity 

Through the above two numerical tests, we clearly 
shown that our method could be able to efficiently and 
exactly deal with the dipole-dipole interaction. In the 
following of this section, dipole-dipole interaction in pho- 
tonic crystal cavity (PCs) has been numerically studied. 
The sketch of the photonic crystal cavity is displayed 
in Fig. |U It consists of a thin dielectric slab with air 
holes arranged as triangular lattice. The lattice constant 
is a, the radii of the air hole is r = 0.3a and the slab 
height is d — 0.6a. The refractive index of the slab is 
3.4. There is a defect hole in the center with refractive 
index n^e/ = 2.4. Its radii is the same as the air hole. 
The two dipoles are located at point A and B on the 
center plane of the slab and the two dipole moments are 
parallel to the x axis. The origin of the axes is set at 
the center of the PCs cavity. For convenience, the spe- 
cial points are also drawn in Fig. 0] 'Atom' A located 
at (0, -11/15, 0)a, 'Atom' B located (0,R- Ho/15,0), 
where R represents the separation and varies from a/ 15 
to 50a/ 15. The two transition dipole moments are paral- 
lel to the x axis. Thanks to the scaling law, the cooper- 
ative decay parameters Tij{u), dipole-dipole interaction 
potential Ay (u>) and frequency U) in unit of 2-nc/a arc 
dimensionless. 

The frequency dependent characters for cooperative 
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FIG. 5: (Color online). The frequency dependent characters 
for cooperative decay parameters and dipole-dipole interac- 
tion potential, (a) The cooperative decay parameters Tij(ui) 
and (b) the dipole-dipole interaction potential Aj 3 -(w) versus 
ui for atom A located at (0, —11/15, 0)a, atom B located at 
(0,11/15, 0)a in the photonic crystal nanocavity. Tij(uj) and 



Ajj(cj) are in unit of oca (a. 



fcTiy/ctp/ncI a). The to is 



in unit of 2nc/a. The two dipole moments are parallel and 
along the x axis. The insets show the behavior around the 
defect frequency (uj c — 0.3133(2-7rc/a)) of the cavity. The two 
transition dipoles are parallel and along the x axis. 



decay parameters Tij(uj) and dipole-dipole interaction 
potential Ay(w) are presented in Fig. [SJa) and Fig. 
djb) for atom A located at (0, —11/15, 0)a and atom 
B located at (0, 11/15, 0)a. The insets show the be- 
haviors for frequency around the resonance frequency 
(cj c = 0.3133(27rc/a)). From Fig. [5ja), we clearly see 
that Tij(ui) oscillates remarkably when u> is far away 
from u c , which shows the powerful modulation ability 
of photonic crystal for the electromagnetic eigen mode 
and the density of photonic states. For uj = uj c (the 
inset), Tij(uj c ) m 144Fo for relatively large atomic sep- 
aration (R 0.46Ao), where r is the largest Tij(oJo) 
for two dipoles located in vacuum for the atomic separa- 
tion R = 0. This large is attributed to the enhance- 
ment of photon emission and reabsorbing rates, which 
can be roughly characterized by the ratio of Q/V. Aimed 
for quantum computation and quantum information pro- 
cessing HI, 59[ where large Fy or large Ay is needed, we 
can increase the quality factor Q through improving the 
nanocavity design. For ui ~ lo c (the inset), Tij(ui) varies 
sharply and is sensitive to the frequency. This implies 
that tuning uj c of the cavity or uii of the dipole can both 
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help to control . Figure EJb) shows similar proper- 
ties for the dipole-dipole interaction potential. The inset 
shows that repulsive or attractive potential can be got 
around the cavity frequency. 



100 
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FIG. 6: (Color online). The position dependent character- 
istics for cooperative decay parameters and dipole-dipole in- 
teraction potential, (a) The absolute maximum cooperative 
decay parameters T*j = and (b) the absolute maximum 
dipole-dipole interaction potential A*j = |A;j versus the in- 
teratomic separation R for atom A located at (0, —11/15, 0)a 
and atom B located at (0, R — lla/15, 0) in the PCs cavity. 
The maximum are calculated for uj around uj c . r*,- and A*j 
are in unit of a.ij (ctij = y / o7y / aJ27rc/o) and ui is in unit of 
2-Kc/a. The two transition dipoles are parallel and along the 
x axis. 



The position dependent characters have been inves- 
tigated. Figure HJa) and [6jb) show the absolute max- 
imum value for lo around lo c of cooperative decay pa- 
rameters r*j = (r^'l and dipole-dipole interaction po- 
tential A*j = | Ay | as a function of their separation R 
(atom A located at (0, —lla/15, 0) and atom B located 
at (0,i? — lla/15, 0)) are show where the maximum are 
calculated for u) around lu c . F*j and A*^ are in unit of 
Oij (ctij = v /757 v /o7J27rc/a). The ui is in unit of 2nc/a. 
The maximum value r*. and A|- varies in a similar way 
except for the amplitude. The above phenomena can be 
understood as follows: r*- depends much on the electric 
field of the defect mode. Equation (fT2|) shows that T*j 
changes in the same way as the electric field of the cavity 
mode at the location of atom B. The cavity mode along 
the y axis looks much the same as Fig. Uta) and we did 
not show it here. Further, we clearly see that T*j = T* u 
and A*j = AJ have been realized for atom A and atom 
B located at positions with same electric field strength 
along the x axis. 



V. CONCLUSION 

In summary, we have proposed a new efficient and rig- 
orous numerical method to investigate dipole-dipole in- 
teraction in photonic crystal nanocavity. We calculate 
the collective and individual radiation rates of classical 
dipoles by directly solving Maxwell's equations in real 
space with a free-space boundary condition. By using 
the result of two dipoles radiation rate minus the sum of 
the two individual radiation rates, we get the cooperative 
decay parameters and dipole-dipole interaction potential. 
Through the self-consistent procedure (Eq. ([6])), exact 
non-Markov results can be got and both real and virtual 
photon effects have been taken into account. Further, it 
can be applied for dipoles with different transition fre- 
quencies in both weak and strong coupling regimes. Our 
investigation suggests that this method works well for 
dipoles located in media with arbitrary shape and may 
be generalized to calculate many dipoles interaction. Nu- 
merical validation has been made in vacuum and planar 
nanocavity. Both results agree very well with the ana- 
lytic. 

Applying this method to a simple photonic crys- 
tal nanocavity, it is found that the cooperative decay 
parameters and the dipole-dipole interaction potential 
strongly depend on the atomic position, the transition 
frequency, quality factor and the cavity frequency. For 
two dipoles arranged at positions with the same local 
coupling strength, the cooperative decay parameters is 
equal to the local coupling strength. Large cooperative 
decay parameters is achieved at the resonance frequency. 
Dipole-dipole interaction potential changes continuously 
from attractive to repulsive case for transition frequency 
varying in a domain of the cavity linewidth around the 
resonance frequency. Larger value and sharper change of 
cooperative parameters and dipole-dipole interaction can 
be obtained for higher quality factor. Owing to the Pho- 
tonic crystal nanostructure is one of the most promising 
platform. It offers many advantages, including high qual- 
ity factor, small mode volume, integrability with waveg- 
uide and fixed position of dipoles. In addition, many 
ingenious schemes for static and ultra fast dynamic con- 
trol of the quality factor, resonance frequency have been 
proposed. Based on these techniques for photonic crys- 
tal nanocavity and methods for tuning the transition fre- 
quency of quantum dots, our results provide some ma- 
nipulative approaches for dipole-dipole interaction with 
potential application in various fields such as quantum 
computation and quantum information processing based 
on solid state nanocavity and quantum dot system. 
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